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Abstract 

We present new evaluation of the critical exponents of 0(ra)-symmetric 0^ 
theory from the field theoretical renormalization group, based on the new 
algorithm for summing divergent series. The central values practically coin- 
cide with those by Le Guillou and Zinn- Justin (1980) but their uncertainty 
is essentially smaller. 
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1. Introduction 

The present paper has an aim to give new evaluation of the critical exponents of 0(77,)- 
symmetric (p^ theory from the field theoretical renormalization group (RG) [IJ, based on 
the new algorithm for summing divergent series [2j. 

According to the formalism of the field-theoretical RG, one should calculate three func- 
tions (3{g), f]{g), V2{g) entering the Callan-Symanzik equation, find a non-trivial root g* 
of equation P{g) = (determining the fixed point of the RG equations), and then the 
critical exponents t] and u, as well as the exponent u of correction to scaling, are given by 
expressions 

V = vi9*), i^-' = 2-r]ig*)+r]2ig*), ^ = /?'(/) • (1) 

The RG functions are given by factorially divergent series in powers of the coupling constant 
g and to calculate them one need a method for summing divergent series. The examples 
are Pade-Borel [3] or conformal Borel [1] techniques. 

Our initial information is the same as in the paper [S], i.e. the first 7 expansion co- 
efficients of the RG functions P{g), vid): V^ig) [31 E] and their large order behavior ^ 
established in the framework of the Lipatov method [7j. The main difference from the 
preceding papers consists in the fact that explicit interpolation of the coefficient function is 
made from the very beginning: the low order expansion coefficients are smoothly interpo- 
lated with their large order asymptotics and unknown intermediate coefficients are found in 
a certain approximation. Considering these coefficients as exact, one can sum the divergent 
series with (in principle) arbitrary precision. The only uncertainty of the algorithm is re- 
lated with ambiguity of interpolation, which has a clear physical sense and originates from 
incompleteness of initial information. As a consequence, the relation of the summation re- 
sults with the assumed behavior of the coefficient functions can be constructively analyzed 
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and estimation of uncertainty becomes completely transparent. For technical reasons, such 
procedure was impossible in the conventional algorithms due to catastrophic increase of 
errors in the course of the series resummation, which made interpolation to be useless. A 
crucial point is stability of our algorithm with respect to smooth errors, involving ambiguity 
of interpolation. 



2. Summation procedure 

Our summation procedure P] is based on the fact that the divergent series 

oo 

w{g) = E w^i-gf (2) 

N=No 

whose coefficients have asymptotic behavior W^^ = ca'^r{N + b), after the Borel transfor- 
mation 

W{g) = J dxe~^x'^^-'B{gx) , B{z) = B^{-zf , = (3) 



N=No 



(where 6o is an arbitrary parameter) and conformal mapping^ z = u/ {1 — u)a reduces to a 
convergent power series in u with coefficients 

N n 

Uo = Bo, UM=j:-;^i-^fC§Zl (iV>l), (4) 



K=l 



whose asymptotics at ^ oo 

Un = UooN--' , = (5) 

a°r(a)r(oo + a) 

is related with the strong coupling asypmptotic behavior of the function W{g), 

W{g) = W^g'^ (g^oo). (6) 

The coefficients Un for iV < 40 are calculated straightforwardly by Eq.(4) and then are 
continued according to power law (5) in order to avoid the catastrophic increase in errors 
[2]. Thus, all the coefficients of the convergent series are known and this series can be 
summed with (in principle) arbitrary accuracy. This completely removes the problem of 
the dependence of the results on variation in the summation procedure, which is the main 



^ This conformal mapping is different from that used in [H [5] . Its advantage consists in the more slow 
growth of random errors in the coefficients Un of resummed series (4) and " super-stability" of the algorithm 
with respect to smooth errors [2]. 



2 



disadvantage of the commonly accepted methods. The interpolation is performed for the 
reduced coefficient function 



Wn , Ai A2 Ak 

= ttt;:! = 1 + + ... + • • • + ... .-..^ + • • • (7) 



W^' N-N {N-Ny (N-N) 



by truncating the series and choosing the coefficients Ak from the correspondence with 
the known values of the coefficients Wl^, Wlq+i, • • •, Wl- The Lipatov asymptotics W^^ 
is taken in the optimal form W^'' = ca^ N''~^^'^T{N + 1/2) [2j[§, and the parameter N 
is used to analyze uncertainty in the results. The Lq value sometimes does not coincide 
with A^o appearing in Eq.2. Indeed, the coefficient function Wn continued to the complex 
plane has a singularity at the point N = a, where a is the exponent of the strong-coupling 
asymptotics (6) [2J. If the exponent a is larger than A'^o, the interpolation with the use of 
all the coefficients is inapplicable: it is necessary to set 

W{g) = WNog""' + ■■■ + W^.g""^ + W{g) , N, = [a] , (8) 

produce summation of the series for W{g), and add the contribution from the separated 
terms; thus, the value [a] + 1 ([. . .] is the integer part of a number) is taken for Lq. Analysis 
of the two-dimensional case [5] shows that a is larger than A^q for almost all the functions. 

One can see that realization of this programe includes (at the intermediate stage) deter- 
mination of the strong coupling asymptotics of the function W{g). A large accuracy of this 
asymptotics is not necessary for summation in the region g ^ 1 and its rough estimation 
is sufficient. 

Following the tradition, we sum the series not only for the functions P{g), rj{g), rj2{g), 
but also for the functions v'^^g) = 2 + r]2{g) — vid) 7~^(5') = 1 ~ '72(fi')/(2 — vid)) 
order to verify self-consistency of the results. A set of possible interpolations was restricted 
by two natural requirements [8j: (a) the interpolation curve comes smoothly through the 
known points and does not have essential kinks for non-integer A^; (b) large A^ asymptotic 
behavior is reached sufficiently quickly, and nonmonotonities at large A^ are on the same 
scale, as a relative difference of the last known coefficient from the Lipatov asymptotics. 



3. The polymer case {n = 0) 

Initial information is given by expansions [3], E] 

p(g) = -g^g^ - 0.4398148149^^ + 0.3899226895^^ - 0.4473160967/ + 0.63385550/ 

-1.034928^^ + . . . + ca^T{N + h)g^ + . . . , 
T]{g) = (1/108)/ + 0.0007713750/ + 0.0015898706/ - 0.0006606149^^ (9) 

^We tried another parametrizations of the form W^'^ — ca^N^T{N + b — b) but the results were 
practicaUy the same, if the same principle was used for restriction of the set of interpolations. 
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+0.0014103421/ - 0.001901867^^ + . . . + c'a^T{N + h')g^ + . . . , 
r^^[g) = -{1/A)g + (1/16)^2 - 0.0357672729^^ + 0.0343748465^ - 0.0408958349/ 
+0.0597050472/ - 0.09928487c/^ + . . . + c"a^T{N + b)g^ + ... , 
with the parameters [6j 

= 0.16624600, 6 = 6' + 1 = 4 , c = 0.085489 , c' = 0.0028836 , c" = 0.010107. 

(10) 

Below we discuss some technical details of the summation procedure. 

Function (3{g). All the interpolations with Lq = 1 are unsatisfactory: the interpolation 
curves rapidly achieving their asymptotic behavior exhibit a sharp kink in the interval 
1 < < 2, indicating a singularity in this interval. Estimation of the strong-coupling 
asymptotics yields a ~ 1, confirming the singularity at ^ 1 and indicating that the 
choice Lq = 2 should be made. In this case, the interpolation curves with A^ < —0.9 
exhibit significant nonmonotonicity at large A^, and the curves with A^ > 1.1 have a kink 
in the interval 2 < A^ < 3 (see Fig.l). Thus, the "natural" interpolations correspond to the 
interval —0.9 < N < 1.1. The summation results are shown in the inset in Fig.l, which 
indicate that 

/ = 1.420 + 1.426, cj = 0.784 ^ 0.795 . (11) 

The g* value is in agreement with the results of early works {g* = 1.421 ± 0.004 \3\, g* = 
1.421+0.008 |1]) and in a certain conflict with the more recent evaluation g* = 1.413+0.006 

Function i]{g). According to Eq. (3), the expansion for ri{g) begins with g'^. We fail 
to obtain satisfactory interpolations with Lq = 2: the curves rapidly approaching large A^ 
asymptotic behavior exhibit a kink in the interval 2 < A^ < 3, indicating that the exponent 
a lies in the same interval. Indeed, the estimate of strong-coupling behavior gives a ^ 2 
and suggests the choice Lq = 3. In this case, the satisfactory interpolation curves (see 
Fig.2,a) exist only for 1.5 < A^ < 2.2. They could be considered unsatisfactory due to 
a kink for 3 < A^ < 4; however, the curves of such a shape provide the exact t] value 
in the two-dimensional case [8]. In our opinion, such interpolations are allowable because 
the amplitude of oscillations of the coefficient function is on the order of the amplitude 
of oscillations of the known coefficients. The summation results are shown in the inset in 
Fig.2,a and give r] = 0.0269 0.0275. 

Functions ri2{g), v^^{g) '-/^^{g). Rough estimates of the strong coupling behavior for 
772(5') and u'^lg) gives a ~ 2 but in general the results are not self-consistent and violate 
the relations between RG functions. Analysis of the 2D case ^ shows that it is related with 
specific features of 77(5): due to small expansion coefficients, this function is small for g^lO, 
but grows rapidly for large g. As a result, asymptotic behavior of i^'^ig) and i]2{g) contains 
a mixture of g and g^ terms, which is difficult to analyze numerically. Therefore, summation 
of the series for 772(5) and i>~^{g) is performed at Lq = 3|1 in order to take into account a 

^ Summation at Lq = 2 gives practically the same results but with lesser uncertainty. 
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Figure 1: Interpolation curves for the expansion coefficients of P{g) and summation results 
for g* and lo. 
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Figure 2: Interpolation curves for the expansion coefficients of functions ri{g) (a), rj2{g) (b), 
i^~^{g) (c) and 7~^(5') (d). The insets show the summation results at g = g*. 
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possible singularity at ~ 2, while the series for 'y~^{g) is summed at Lq = 1, but without 
the restriction of kinks for noninteger N (due to relation j'^lg) = 1 + r]2{g)/{2 — i]{g)) 
its coefficient function is expected to be regular for > 1, but containing a smeared 
singularity at ~ 1). Figures 2,b-d show the allowable interpolations and summation 
results. The latter are presented in Table 1 and compared with results by other authors 



T a b I e 1. 

Critical exponents for the polymer case (n = 0) from the field theory 





BNM 


LG-ZJ ^ 


G ZJ 


Kl dU] 


j-Ki [m 


Present work 


7 


1.161(3) 


1.1615(20) 


1.1596(20) 


1.161 


1.1604(8) 


1.1615(4) 


u 


0.588(1) 


0.5880(15) 


0.5882(11) 


0.5883 


0.5881(8) 


0.5886(3) 


V 


0.026(14) 


0.027(4) 


0.0284(25) 


0.0311(10) 


0.0285(6) 


0.0272(3) 


V2 


-0.274(10) 


-0.2745(35) 








-0.2746(7) 


UJ 


0.794(6) 


0.800(40) 


0.812(16) 


0.810 


0.803(3) 


0.790(6) 


g* 


1.421(4) 


1.421(8) 


1.413(6) 






1.423(3) 



The results for u obtained by summation of different series (in view of relations 7 = 
i>{2 — ■)]), = 2 + ri2 — r],i'= (1 — 7)/?72) are presented in Table 2. The fourth estimate is 
rather inaccurate and will be ignored, while the first three estimates practically coincide; 
the relative shift of the central values for them can be considered as the scale of the non- 
controllable systematic error, 

6syst^ 0.0002, (12) 

which appears because the "natural" interpolations for different interdependent functions 
are not completely consistent. For the two-dimensional case [8], this effect is the main 
source of the error: a similar estimate gives Sgyst ~ 0.05, which is larger than the natural 
summmation error for most functions. At the present case 6syst is rather small. 

T ab I e 2. 
Different estmates for the exponent u 
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Series 


Interval for v 


Central value 




0.5883^0.5889 


0.5886 




5885 — 5891 


5888 


V2{9), vig) 


0.5884^0.5892 


0.5888 




0.5848^0.5916 


0.5882 



4. The Ising universality class (n = 1) 

Initial information is given by expansions [3], [5] 

p(g) = -g + g^ - 0.4224965707^^ + 0.3510695978^^ - 0.3765268283/ + 0.49554751/ 

-0.749689^^ + . . . + ca^T{N + b)g^ + . . . , 
ri(g) = (8/729)^2 + 0.0009142223/ + 0.0017962229/ - 0.0006536980/ (13) 

+0.0013878101/ - 0.001697694^^ + . . . + c'a^T{N + h')g^ + . . . , 
r^^[g) = -{l/3)g + (2/27)^2 - 0.0443102531/ + 0.0395195688/ - 0.0444003474/ 
+0.0603634414/ - 0.09324948^^ + . . . + c"a^T{N + b)g^ + . . . , 
with the parameters [6] 

a = 0.14777422, 6 = 6' + 1 = 4.5 , c = 0.039962 , c' = 0.0017972 , c" = 0.0062991 . 

(14) 

The situation is qualitatively analogous to the previous case, and we use the same values 
for the parameter Lq, i.e. Lq = 1 for ^'^{g), Lq = 2 for l3{g), Lq = 3 for other functions. 
Admissible interpolations correspond to the intervals —1.0 < N < 1.2 for j3{g), 1.6 < N < 
2.2 for r/(/, -0.5 < iV < 2.1 for ly'^g), -6.2 < iV < 2.6 for r]2{g), -1.1 < N < 0.95 for 
"y~^{g), and the appearance of the interpolation curves is visually close to that for a case 
n = (Figs. 1,2). The results are presented in Table 3. 

Tables. 

Critical exponents for the Ising case (n = 1) from the field theory 
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BNM [3] 


LG-ZJ [4J 


G-ZJ 


Kl [lU] 


j-Ki [m 


Present work 


7 


1.241(4) 


1.2405(15) 


1.2396(13) 


1.241 


1.2403(8) 


1.2411(6) 


V 


0.630(2) 


0.6300(15) 


0.6304(13) 


0.6305 


0.6303(8) 


0.6306(5) 


V 


0.031(11) 


0.032(3) 


0.0335(25) 


0.0347(10) 


0.0335(6) 


0.0318(3) 


Tlo 
'12 


-0 382f5) 


-0 3825f30) 








-0 3832f8) 




0.788(3) 


0.790(30) 


0.799(11) 


0.805 


0.792(3) 


0.782(5) 


9* 


1.4160(15) 


1.416(5) 


1.411(4) 






1.4185(25) 



5. The XY universality class {n = 2) 

The detailed discussion of this case is given in the paper [9]. For completeness, we 
present here the final results of this study (Table 4). 

T a b I e 4- 

Critical exponents for the XY (or helium) case (n = 2) from the field theory 





BNM [3] 


LG-ZJ [4J 


G-ZJ [5] 


Kl do] 


J-Kl [llj 


Present work 


7 


1.316(9) 


1.3160(25) 


1.3169(20) 


1.318 


1.3164(8) 


1.3172(8) 


u 


0.669(3) 


0.6695(20) 


0.6703(15) 


0.6710 


0.6704(7) 


0.6700(6) 


V 


0.032(15) 


0.033(4) 


0.0354(25) 


0.0356(10) 


0.0349(8) 


0.0334(2) 


V2 


-0.474(8) 


-0.4740(25) 








-0.4746(9) 


LU 


0.780(10) 


0.780(25) 


0.789(11) 


0.800 


0.784(3) 


0.778(4) 


9* 


1.406(5) 


1.406(4) 


1.403(3) 






1.408(2) 
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6. The Heisenberg universality class (n = 3) 
Initial information is given by expansions [Sj [5] 

p(g) = -g + g^ - 0.3832262015^^ + 0.2829466813c/^ - 0.27033330c/^ + 0.3125559/ 

-0.414861c/^ + . . . + ca^T{N + b)g^ + . . . , 
r^{g) = (40/3267)c/2 + 0.0010200000^ + 0.0017919257c/^ - 0.0005040977^^ (15) 
+0.0010883237/ - 0.001111499c/^ + . . . + c'a^T{N + h')g^ + . . . , 
r^^[g) = -(5/11)^ + (10/121)/ - 0.0525519564/ + 0.0399640005/ - 0.0413219917/ 
+0.0490929344/ - 0.06708630^^ + . . . + c"a^T{N + b)g^ + ... , 
with the parameters ^ 

= 0.12090618, 6 = 6' + 1 = 5.5, c = 0.0059609 , c' = 0.0003656 , c" = 0.0012813 . 

(16) 

The same values for Lq, as in previous cases, were used. Admissible interpolations corre- 
spond to the intervals -1.0 < iV < 1.6 for p{g), 1.6 < N < 2.3 for r]{g) OA < N < 2.0 for 
i/~^((y'), —0.6 < N < 2.2 for 7/2(5'), 0.5 < N < 0.95 for 'j~^{g). The results are presented in 
Table 5. 



T a b I e 5. 

Critical exponents for the Heisenberg case {n = 3) from the field theory 





BNM [3] 


LG-ZJ [4J 


G-ZJ 


Kl [TU] 


J-Kl [TT] 


Present work 


7 


1.390(10) 


1.386(4) 


1.3895(50) 


1.390 


1.3882(10) 


1.3876(9) 


u 


0.705(5) 


0.705(3) 


0.7073(35) 


0.7075 


0.7062(7) 


0.7060(7) 


V 


0.031(22) 


0.033(4) 


0.0355(25) 


0.0350(10) 


0.0350(8) 


0.0333(3) 


V2 


-0.550(12) 


-0.5490(35) 








-0.5507(12) 




0.780(20) 


0.780(20) 


0.782(13) 


0.797 


0.783(3) 


0.778(4) 


9* 


1.392(9) 


1.391(4) 


1.390(4) 






1.393(2) 
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7. Discussion 



One can see from Tables 1,3,4,5 that there is a good correspondence between the different 
field theory estimations. A surprisingly good agreement takes place between our results and 
estimation by Le Guillou - Zinn- Justin [4J: the typical difference of the central values is less 
than 0.0010, in spite of rather conservative estimation of errors given in [1]. This coincidence 
is not in any degree incidental: the authors of [1] carried out interpolation of the coefficient 
function in order to predict one or two of unknown expansion coefficients and used them to 
give some kind of the expert prediction, but were induced to allow rather large uncertainty 
of results due to their strong dependence on variation of the summation procedure. On the 
other hand, recent reevaluation in [5J looks somewhat artificial and has a tendency to shift 
the results beyond their natural range; in particular, the shift of g* in comparison with 
[1] is always made in the direction opposite to ours (Tables 1,3,4,5)0- 

A good agreement 

can be seen also with variational perturbation theory p^; it is especially pleasant that 
taking into account the large order perturbation behavior and more elaborated estimation 
of errors [TT] makes the results more close to ours. A small disagreement still remains for 
the exponent r] but it is on the same level as violation of the relation 7 = z/(2 — rj) for the 
central values of [TT]. 

Now let us discuss the correspondence of our results with other information on the 
critical exponents, provided by physical experiment, Monte Carlo simulations (MC) and 
high temperature series (HT) [12]. 

Case n = 3. Overall scattering of the MC and HT results is rather large and in this 
extent they agree with Table 5. There is a tendency to a small disagreement between our 
and the most recent MC resufts (7 = 1.3960(10), z/ = 0.7112(5) [13]) but the latter are 
in the same disagreement with the physical experiments, whose results for 7 are grouped 
around value 1.386 (see Tables 24,25 in [i2\). Analogously, the experimental results for 
the exponent (3 suggest the mean value 0.365 [12] in the good agreement with our estimate 
f3 = 0.3648(4) (following from Table 5) and in the worse agreement with value f3 = 0.3689(3) 
given in [13] . 

Case n = 2. A situation is analogous to the previous case. Overall scattering of the 
MC and HT results is rather large (see Fig.l in ^9j) but the recent results have a tendency 
to contradict Table 4 (7 = 1.3178(2), u = 0.6717(1) r] = 0.0381(2) [H]). Simuftaneously 
they contradict the experiments in liquid helium, i.e. value 



obtained by the measurements of superfiuid density from the second sound velocity [15] . 
and the results 



u = 0.6705 ±0.0006, 



(17) 



a = -0.01285 ±0.00038, 



u = 0.67095(13) [IS] , 



a = -0.01056 ±0.00038, z/ = 0.6702(1) [IT], 



(18) 



^It should be noted that [HIS] and the present paper use the same information for (3{g). 
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a = -0.0127 ±0.0003, 1^ = 0.6709(1) [IE] 

obtained in the satellite measurements of the thermal capacity (the relation a = 2 — du 
was used). 

Case n = 1. In this case, the HT and MC results are numerious (see Tables 3,5 in [12j ) 
and can be summarize as 

7 = 1.2372(5) 

u = 0.6301(4) (19) 
r] = 0.0364(5) 

(see Eq.3.2 in [12]). One can see from the Table 3 that beautiful consensus was reached 
for the exponent z/; on the other hand, values for 7 and rj in (19) are in the meaningful 
contradiction with Table 3. The experimental results have large uncertainty and cannot 
compete with theoretical predictions. 

Case n = 0. In this case, precise results for the exponent u can be obtained by direct 
study of self-avoiding walks on the lattice; due to simplicity of the algorithm, a good 
statistics can be gathered. The most recent results (z/ = 0.5876(2) \19\, v = 0.5874(2) 
|20j . ly = 0.58758(7) pT]) appear in a slight contradiction with our result in Table 1. This 
contradiction is not very significant and we can avoid it by allowing more wide set of 
interpolations and extending the error bars; however, somewhat "unnatural" interpolation 
curves should be used for it. The results for the exponent 7 are essentially less precise [12] 
and cannot compete with Table 1. 

We can conclude that the general situation is satisfactory but there is a cause for anxiety 
related with the exponents 7 and rj in the Ising case. Disagreement on the scale 0.003 is 
essentially larger than uncertainty of the recent field theoretical estimations (Table 3) and 
that of Monte Carlo results. It is also essential that the latter are obtained by different 
researches and not related with a specific group. At present, the origin of this disagreement 
is unclear and further investigations are necessary. 

This work is partially supported by RFBR (grant 06-02-17541). 
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